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!• The physical probl c/i?: . 

We considered a flow about an oscillating thin airfoil in 
a transonic stream. We assume that the flow field can be 
decomposed into a mean flow plus a periodic perturbation. 

Letting cp(x,y) denote the velocity potential associated with 
the latter we have 

(1) ([K - (y-l)Uo]cPx3x '^yy “ (2iw/e)cp^ 

+ [u)Ve - iu)(y-l)u^^]cp = 0, 

where K = {1-M )/M e, M is the free stream Mach number, y is 
the ratio of specific heats, e = (6/M)^^^, 5 is a measure of the 
airfoil thickness, UQ(x,y) is the mean flow velocity, x anc, y 
have been suitably scaled, and cp is the velocity potential of 
the flow perturbation. On the surface of the airfoil the usual 
Neurnman conditions arc imposed. For the problem considered here, 
this boundary condition takes the form 

{2) cp = (f,) + iw(f,) for y = +0, x on the wing 

y X X. """ 

where f^(x) and f_(x) described the spatial behavior of the 
upper and lower surfaces of the wing, respectively. 

Since the mean flow is transonic, (1) is a partial differen- 
tial equation of the mixed type. Furthermore, for large enough 
values of the frequency, (1) is an indefinite equation. If in 
the far field the flow is uniform and subsonic, (1) reduces to 

(3) (l-M^)cp^^ - 2ituMcp^ + cp^ + (u^cp = 0, 

which is a convected Helmholtz equation. 


2 


Of crucial importance is the correct radiation at infinity, 
M. Gunzburger* has shown that this is given by 

.2 


(4) 




where R = (x^ -h (1 - is large. 

To complete the specification of the problem we need to 
impose a condition across the wahe. Following standard practice 
we assume that the wake is a line starting at the trailing edge 
of the airfoil. Continuity of pressure and normal velocity are 
required^ and assuming the airfoil lies on the line y = 0, these 
become 


(5) 


where 


[cp^ t iu)cp] = O 


[«Py] = 0 


y = 0, X in the wake region 


(6) [%] = |(x,0-!-) - 5(x,0~). 

The problems of numerical approximation are tv\?o fold.** 

The first and more sei-ious one is that for a range of tu of 
physical interest, the matrix problem resulting from a discriti- 
zation of the partial differential equations is indefinite and 
non-symmetric . This has made the use of standard iterative 
techniques for solving this sparse matrix system impossible to 
use. The second difficulty is that different descritizations 

^See Final Report on NASA contract NSG-1366, 

H, Weatherill, F. E. Ehlers, and J. D. Sebastian, NASA CR-25 99, 
1975. 


had to bo used in the hyperbolic region and the elliptic region. 
The new method outlined in the next section obvintes both of 
these difficulties. 

2 . The least squairos finite element metho d . 

As a first step we reduce the governing system (1), (2), (3) 
and (4) into an equivalent first order system 

([K - (V + ■ (2iui/e)u + (uiVe - io)(v-l)Ug^]cp = 0 

- u = 0 



V = f^ -h iiuf on the wing 
Au -l- iuiAcp = 0 

on the wahe 

Av = 0 

. 2 

^ u + ^ V ‘f* (1 - ^)cp = 0 in the far field. 

1-M 

Gunzburger's radiation condition (3) will be imposed on a rectangle 
enclosing the wing, and whose edges are "far removed" from the 
wing. The least squares method consists of minimizing, over a 
finite element space of vector valued fvinction, a fv;nctional 
which is obtained by taxing the square of the absolute value of 
each of the above terms and integrating over the appropriate . 
region. 
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This yields a linear system of algebraic equations whose 
coefficient matrix is banded . sparse . Hermitian and positive 
definite . Therefore iterative methods , which take full advantage 
of the sparcity, can be used. Furthermore, the method uses the 
same discritization technique in elliptic and hyperbolic regions . 

3 . Implementation of method . 

Two computer programs have been written, both using linear 
basis functions over triangles for the finite element space. The 
first program uses a banded Gaussian elimination solver to solve 
the matrix problem, while the second uses an iterative technique, 
namely SOR, Both programs are considered to bp working, but 
certainly not optimized. The only results obtained are for an 
oscillating flat plate (not a true transonic problem) although the 
programs have been written to admit general input. The major 
accomplishment of these programs is their success for values of 
(u, using SOR, well above the critical value at which other 
approaches failed. 

The following projects would be of immediate interest: 
a) Improving the accuracy of the codes by using different grids. 
Theoretical and numerical experiments on model problems have 
shown that a triangularization of the type (i) below yields 
answers on order of magnitude more accurate than that using 
type (ii) (which is presently being used in the program) . It 
is a relatively triv’;lal matter to convert the program to type 
(ii) triangles. 
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(i) (ii) 

b) Inputing some realistic transonic mean flows. The programs 
are built to accept such input, so this is also a relatively easy 
matter. 

c) Optimize the codes to achieve maximum efficiency. 

d) Write an extensive report detailing the method and results of 
its implementation. 

Of longer range interest are the solution of 3-D problems ^ 

and of problems with complicated geometries . In both cases the 
difficulty of the task is alleviated by the least squares method, 
in the first case because of the ability to use iterative methods, 
and in the second by the flexibility of finite element methods. 



